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Oh Keywords. Fast Diffusion, Total Variation Flow, asymptotic, extinction profile, convergence 

■^ rates, extinction time. 

^' 

B 

(N 
> 

cn 
in 

(N 

l> 

O 



X 

(a) Departamento de Matematicas, Universidad Autonoma de Madrid, Campus de Canto- 
bianco, 28049 Madrid, Spain. 

E-mail: niatteo.bonforte@uam.es. Web-page: 'http : //www . uam . es/matteo . bonf orte' 

(b) Department of Mathematics, The University of Texas at Austin, 1 University Station 
C1200, Austin, TX 78712-1082, USA. 



E-mail: figalli@math.utexas.edu. Web-page: http://www.ma.utexas.edu/users/figalli/ 



1 Introduction 

The Total Variation Flow (TVF) and the Sign-Fast Diffusion Equation (SFDE) are two de- 
generate parabolic equations arising respectively as the limit of the (parabolic) p-Laplacian as 
p ^ 1+ and of the fast diffusion equations as ?7i — t- O"*". In one spatial dimension these two 
equations are strictly related, and our goal here is to describe their dynamics. 

Let us first introduce the SFDE: consider the Fast Diffusion Equation (FDE) 

dtv = A{v"'), 0<m<l (1.1) 

(by definition, v"^ = \v\'^~^v). By letting m — )• 0"*" one gets the SFDE 

dtv = A( sign(f)) . 

To study the evolution of this equation in one dimension, we will exploit its relation with the 
TVF (also called 1-Laplacian, as it corresponds to the limit of the p-Laplacian when p — )• 1^"): 
at least formally, if v solves the SFDE, then u{x) := L v solves the TVF 

Of course this is purely formal, and we will need to justify it, see Section [3J The correspon- 
dence between solutions of the p-Laplacian type equations and solutions to fast diffusion type 
equations have been used since a long time, cf. [23] and more recently in [6], in the study of 
equations related to the fronts represented in image contour enhancement. In several dimen- 
sion the correspondence between solutions of the p-Laplacian and of the FDE is less explicit 
and holds only for radial solutions, cf. [2^. Here our strategy is first to analyze the dynamic 
of the TVF in one dimension, and then use this to recover the behaviour of solutions of the 
SFDE. 

The literature on the TVF is quite rich, and we suggest to the interested reader the monograph 
[5] as source of references for the existence, uniqueness, basic regularity and different concepts 
of solutions and their relations, together with estimates on the extinction time (see also the 
review paper |16j and references therein). The TVF has some interest in applications to noise 
reduction, cf. [23l HI E]. The asymptotic behaviour of the TVF is still an open problem in 
many aspects, even if partial results have appeared in [H El [SJ [71 [HI [9l [TTj. However, at least 
in the simpler case of one spatial dimension, the results contained in the present paper exhibit 
an almost explicit dynamic and a sharp asymptotic behaviour. After the writing of this paper 
was essentially completed, we learned of a related work J19j . 

Plan of the paper. 

• In section [2] we analyze the dynamic of the one-dimensional TVF. As a first step in this 
direction, we study the time discretized case, for which we find the explicit dynamics for "local 



step functions", see Subsection 2.3 and 2.4 In Subsection 2.5 we pass to the continuous time 



case, and we find the explicit evolution under the TVF for a generic "local step functions". 



Then, in Subsection 2.6, exploiting the stability of the TVF in L^ spaces and arguing by 



approximation we prove some basic but important properties of solutions to the TVF, such as 



the conservation and contractivity property of the local modulus of continuity (Theorem 2.8), 
and the explicit behaviour around maxima and minima. 



In the case of nonnegative compactly supported initial data, we first prove an explicit formula 
for the loss of mass and extinction time for the associated solution to the TVF (Proposition 



2.10). Next we study the asymptotic behaviour of such solutions: we analyze and classify 



the possible asymptotic profiles (that we shall call more properly "extinction profiles", see 



Theorem 2.11), and we characterize the asymptotic behaviour of solution to the TVF near the 



extinction time as a function of the initial datum (see Theorem 2.15). To our knowledge this 



is the first completely explicit asymptotic result for very singular parabolic equations, even if 



it holds in only one spatial dimension. Finally, in Subsection 2.10 we prove the sharpness of 



the rate of convergence provided by Theorem 2.15 



• In Section [3] we dedicate our attention to the SFDE. First we rigorously show that a BV 
function solves the TVF if and only if its distributional derivative solves the SFDE, and then 
we exploit this describe the dynamic of the SFDE. We conclude the paper with a discussion on 
the relation between the SFDE and the Logarithmic Fast Diffusion Equation (LFDE), which 
is another possible limiting equation of the Fast Diffusion Equation as m — >■ 0^ , see Subsection 
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2 The l-dimensional Total Variation Flow 

In this part we deal with the one dimensional TVF. Before introducing the problem, we first 
recall briefly some notation and basic facts about BV functions for convenience of the reader. 

2.1 Notations and basic facts about BV functions in one dimension 

Here we recall some basic facts about one-dimensional BV functions, referring to ^ Section 
3.2] for more details. 

Consider an open connected interval /CM. A function u £ BV{I) if n G Lj^^{I), its 
distributional derivative Du is a (signed) measure, and its total variation \Du\ has finite mass. 

The distributional derivative Du can be decomposed as Du = dxudx + D^u, where dxU = Vu 
is the absolutely continuous part of Du (with respect to the Lebesgue measure), and D^u is 
the singular part. 

By Sobolev inequalities we have the inclusion BV{R) C L°°(R), and BV{I) n L^(/) C LP{I) 
for any 1 < p < oo. 

If u G BV{I), up to redefining the function in a set of measure zero, for every point x £ I it 
always exists the left or right limit of m at a point x, which we denote by 

u{x ) := lim u{y). 

Moreover, the limits above are equal up to a countable number of points. We will always 
assume to work with a "good representative" , so that the above property always holds (see [21 
Theorem 3.28]). 



2.2 The setting 

Let us briefly recall the definition of strong solution to the TVF, in the form we will use it 
throughout this paper. For the moment, we do not specify any boundary condition (so the 
following discussion could be applied to the Cauchy problem in M, as well as the Dirichlet or 
the Neumann problem on an interval). 

A function u G L°°{[0,oc),BV{I)) n Wl;^{[0,oo),L'^{I)) is a strong solution of the TVF if 
there exists z £ L^^^([0, oo), H^^'^(/)), with ||2;||oo < 1, such that 

dtu = dxZ on (0,00) xl, (2.1) 

and 

/ z{t,x)Du{t,x)dt dx= \Du{t,x)\dxdt VT > 0. 

Jo Ji Jo Ji 

Roughly speaking, the above condition says that z = Du/\Du\. We refer to the book [5] for 
a more detailed discussion on the different concepts of solution to the TVF depending on the 
classes of initial data (entropy solutions, mild solutions, semigroup solution), and equivalence 
among them. 

Throughout the paper we will deal with non- negative initial data for the TVF, although many 
properties maybe extended to signed initial data. 

2.3 The analysis of the time-discretized problem. 

It is well known that the strong solution u of the TVF defined above is generated via Crandall- 
Ligget's Theorem, namely it is obtained as the limit of solutions of a time-discretized problem, 
formally given by the implicit Euler scheme 

u{ti+i) - u{ti) _ r. [ Du{ti+i) 



ti+i -ti \\Du{ti+i)\ 

We refer to the book [5] fore a more complete and detailed discussion of these facts. The goal of 
this section is to understand the behaviour of the time-discretized solution both at continuity 



and at discontinuity points (see Propositions 2.2 and 2.3) 



Let us fix a time step /i > 0, set to = 0, tj+i =ti + h = {i + l)h, and define Uih{x) := u{ih, x) 
so that uo{x) = u{0,x). The first step reads: 

"^^"^7"°^"^ = d.z„ (2.2) 

n 

where Zh G L°°{I) satisfies ||^;/i||oo < 1 , z^Duh = \Duh\, and dxZ G L'^{I)- Of course it suffices 
to understand the behavior of Uh starting from uq, as all the other steps will follow then by 
iteration. 

We are going to prove the main properties of the time discretized solution, and to this end is 
useful to recall an equivalent definition for Uh'- 



Uh 



argmin [$;j («)] , where '^^{u) = / \Du\ + -— / |ti — uqI da; . (2.3) 

J I 2n J I 



Indeed by strict convexity of the functional $ft, the niinimizer is unique and is uniquely charac- 



terized by the Euler-Lagrange equation associated to ^h, which is exactly (2.2), that we shall 
rewrite in the form 

Uh = UQ + hdxZh. (2.4) 

Let us observe that the above construction does not need uq to be L'^{I)'- if uq £ BV{I) the 
above scheme still makes sense and provides a function u^ such that Uh — uq £ L'^i^)- 

Next, we remark that since uq , Uh G BV{I) also dxZh G BV{I) C L°°{I), which implies 
that Zh is Lipschitz and is differentiable outside a countable set of points. Define the (at most 
countable) set 

N(zh) := |x G M lim ^h{x + e) - Zh{x) ^^^^ ^^^ ^^.^ ^ 
[ £-s>o e J 

Since dxZh G BV{I), it is continuous outside N{zh) (i.e. dxZh G C^(R \ N{zh))), and we have 
that N{zh) coincides with the set of discontinuity point of u^ — uq . 



Collecting all the information obtained so far, we can say that equation (2.4) is equivalent to 



hdxZh{x) = Uh{x) — Uo{x) for all x G M \ N{zh) 

\zh{x)\ < 1, for ah X G M (2.6) 

Zh.[x) = ±1 , for \Duh\ — a.e. 

The next lemmata will allow us to show the important fact that, on M \ N{zh), Uh is locally 



constant whenever different from uq (see Proposition 2.2). 
Lemma 2.1 The following holds: 

\Duh\ {{dxZh / 0} \ N{zh)) = \Duh\ {{uh / no} \ N{zh)) = . (2.7) 

Proof. The equality 

{dxZh + 0} \ N{zn) = {uh / no} \ N{zh) (2.8) 

easily follows by observing that u^ — uq = hdxZh G C^{U. \ N{zh)). Moreover, since Zh ■ Duh = 
I Duh I , we have 

/ (l-|z/,(x)|)d|Dn/,|=0 
which in particular implies 

/ {l-\zh{x)\)<l\Duh\=^. (2.9) 

J{Uh^UQ}\N(zh) 

Next we notice that since Zh is differentiable on M \ N{zh) and \zh\ < 1, we have {zh = 



±l}n{dxZh y^O}n{R\N{zh)) = 0. Hence, thanks to ^^ we deduce that {l-\zh{x)\) > on 
the set {uh / no} \N{zh). Combining this information with (2.9) we obtain that \Duh\{{uh / 
no} \ N{zh)) = 0, as desired. [] 

We now use the above lemma to analyze the behavior of u^ near continuity points. 



Proposition 2.2 (Behaviour near continuity points) If Uh is different from uq at some 
common continuity point x, then it is constant in an open neighborhood of x. 

Proof. Let x £ I he a continuity point both for uq and u^, and assume that Uh{x) > uo{x) 



(the case Uh{x) < uo{x) being analogous). Then by (2.6) we deduce that dxZh is continuous 



and strictly positive in an open neighborhood I{x) of x, which together with (2.7) implies 



\Duh\{lix)) =0. 



Hence Uh is constant on I{x). Q 



We now show that if uq G BV{I) has some discontinuity jump, then Uh can only have jumps 
at such points, and moreover the size of such jumps cannot increase. 

Lemma 2.3 (Behaviour at discontinuity points) Let uq G BV{I). Then, the following 
inequalities hold for any x G I: 



if Uh{x ) < n/j(x+) then 

if Uh{x^) < Uh{x^) then 

Moreover, 



uo{x ) < Uh{x ) < Uh{x-^) < uoix+) 
uoix^) < Uhix+) < Uh{x^) < uoix^) ■ 



Uh{x ) < u/j(x~^) implies Zh{x) = 1 

Uh{x^) > Uh{x^) implies Zh{x) = — 1 . 

Proof. Let x G / be a discontinuity point for Uh- Then 

Duh{x) = {uh{x'^) - Uh{x~))5x , 



(2.10) 



(2.11) 



(2.12) 



where 6x is the Dirac delta at x G /. We first prove (2.11): recalling that Zh ■ Duh = LDu/i , 



if Uh{x ) < Uh{x'^) then by (2.12) we get Zh{x) = 1. Analogously Uh{x ) > Uh{x^) implies 
Zhix) = -1. 



Let us now show (2.10): assume first that Uh{x ) < Uh{x^). Since Zh{x) = 1, x is a maximum 
point for z^, thus dxZh{x~) > and dxZh{x^) < 0. Using (|2.6[), this implies 



nfa(x+) --u.o(x+) 
h 



dxZhix^) < < dxZhix ) 



Uh{x ) -uo{x ) 



h 



which combined with our assumption Uh{x ) < Uh{x^) gives 

""0(2;") < Uh{x^) < Uhix"^) < uo{x^). 

The case Uh{x~) > Uh{x'^) is analogous. □ 
As an immediate corollary we get: 

Corollary 2.4 (Local continuity) Let x G /. Ifuo is continuous atx, thenuh is continuous 
at X. 



This result shows that (at least at the discrete level) the TVF cannot create new discontinu- 
ities, and that continuous initial data produce continuous solutions (this is actually what we 
will prove in Theorem 2.8). Observe that this is a local property which does not depend on 
the boundary conditions (that we have not specified yet). 

We conclude this section with the following estimates on the local loss of mass. 
Lemma 2.5 (Local L^— estimates) The following estimates hold for any interval {a,b) C I: 



Uh{x)dx— / uo{x)dx 



< 2h. 



(2.13) 



Proof. By equation (2.6) we have 

rb 

Zhih) - Zh{a) = 



dxZh{x) dx 



1 



Uh{x)dx— / uo{x) dx 



(2.14) 



Hence (2.13) follows from the bound H^hHoo < 1- D 



2.4 The dynamics of local step functions I. The time discretized case 

In this section we use the time discretization scheme to study the dynamics for initial data uq 
which coincide with a step function on some open interval /. 

Let us point out that, if uq is exactly a step function, then one can give an explicit formula 



for its evolution (see Section 2.4.2 and p^) by simply checking that it satisfies the equation. 
However, by studying the "time-discretized" evolution (and then letting /i — )• 0), one can see 
in a much more natural way the "locality" in the dynamic of the TVF. Moreover, our method 
shows how to deal with functions which do not belong to L?{I)). Finally, our description 
give a good insight of the analysis of the discretized PDF, which may be useful for numerical 
purposes. 

We would like to notice that it is important for the sequel that the time step /i > is 
sufficiently small with respect to the size of the jumps of the step function we are considering, 
as otherwise the discretized dynamics becomes more involved, as we shall show with an example 
at the end of this section. 



2.4.1 Local evolution of a single step. 

To give an insight on the way the discretized evolution behaves, we consider the case of maxi- 
mum steps, whose behavior is made clear in Figure [T]. Let us fix an interval I = /i U /2 U I3, 
and assume that uq = aixi + 0^2X2 + osXs on h with 012 > max{ai,a3} and Xk = Xh i^ 
the characteristic function of the open interval Ik = {xk-i,Xk). (Observe that we make no 
assumptions on uq outside /.) Fix h > small (the smallness to be fixed), and consider the 
function Uh- We have: 



(a) Ufi is constant on any interval Ik- This follows easily from by Proposition 2.2 and Corollary 



(b) // h is small enough, then u^ jumps at the points xi and X2- Indeed, if by contradiction 
xi is a continuity point for Uh, Uh would equal to a constant a on /i U I2, with ai < a < a2 




Figure 1: Dynamics of a maximum step. This figure shows the dynamic only inside the interval [xo,X3] 



(see Lemma 2.3). However, by Lemma 2.5 we know that 



and 



which is impossible if 



\ai — a\ |/i 



\oi2 — a| Ih 



[uhix) -^0(2;)] dx 



h 



[uh{x) -uo{x)] dx 



< 2h 



< 2h 



In particular, Uh jumps both at xi and X2 if the "simpler" condition 

< h < \ai — 02! min{|/i| , I/2I} 
holds. 



(c) Applying Lemma 2.3 again (see equation (2.11)) we obtain 

Zhixi) = 1 and Zh{x2) = -1 , 



which combined with (2.14) gives 

rx2 



Uh{x)dx— / uo{x) dx 

Xl J X\ 



X2 



-2. 



Hence 



2/1 



Uh = a2,h ■= "2 



on I2 . 



(2.15) 
(2.16) 



(d) Combining all together we obtain that 



2/1 



uo = aixi + 02X2 + asxs on/ implies w/i = ai.hXi + ( "2 - 777 ) X2 + a3,ftX3 on/, 

^ -'2 



o:2,h 



(2.17) 



where ai < ai^h ^ Q^2./i) Q^3 ^ Q^3,?i. ^ ct2,/i (see Lemma 2.3). The exact value of aih and 03. ^ 

1 al 

2/1 



depend on the behavior of uq outside /, but using (2.14) we can always estimate them: 

2h 



\ai,h - ail < 



|/i| 



\a3,h - 03] < 



l^sl' 



(2.18) 



(In some explicit cases where one knows that value of z at xi and X3, ai^h and a^^h can be 
explicitly computed using (2.14).) 



Remark 2.6 It is important to observe that the value of 02, /» is independent of the value of 



ai^h and a^^h, but only depends on the fact that ai, a^ < 02, /i- In particular, thanks to (2.18), 
one can iterate the above construction: after i steps we get 



Uih = ai/hXi + a2ihX2 + a3PhX3 on /, a2Ph = "2 



2ih 



holds as long as ai,(^-i)/i,a3,(^-i)/i < a2/hi which for instance is the case (by iterating the 
estimate (|2.18[)) if 



£/i < |ai — 02! min{|/i|, II2I} and £/i < |a2 — 03] min{|l2|, II3I}. 

In particular observe that in this analysis we never used that /i and I3 are bounded intervals, 
so the above formulas also holds when xq = — (X) and X3 = +00. 



2.4.2 Evolution of a general step function 

The above analysis can be easily extended to the general A^-step function: assume that 

Af+l 
uo = 



^ akXk on / 



fc=o 



where a^ G M for k = 0, . . . , A^ + 1, and Xk = Xh is the characteristic function of the open 
interval I^ = ixk-i,Xk) (also the values xq = —00 and xat+i = +00 are allowed). Then, if 



< £h < min < \aj — aj+i| min{|/j|, |/j+i|} \ 



(2.19) 



the discrete solution after i steps is given by 



N+l 
Uih = ^ OLk/hXk on /, 
fc=o 



where we are able to explicitly get the values of ak/h for k = 1,...N, see Remark 2.6, and 
some information on ao,ffc and UN+i/k'- for k = 1, . . .N 



Oiklh 



Q-k, 






2th 


Oik - 


\Ik\ 


ak + 


2£h 

141 



if ak-i <ak < ak+i or if Uk+i < ak < Ofc-i 
if ak > max{afc_i , Ofc+i} 
if Ofc < min{afc_i, a^+i} 



(2.20) 



[ < ao,(£-i)/i ' 



OAT+l.ft 



> aAf+i,(£-i)/i ) 

< aAf+l,(£-l)h ) 



if ao < ai 
if ao > Qi 

if a]\f > oat-i-i 

if OiV < OAf+l- 



A concluding remark on the smallness of the time step h. Since we want to describe 
the behaviour of the TVF, we are mainly interested in the hmit /i — t- 0, which means that 
condition (2.19) is always fulfilled. Anyway it is interesting to observe that the dynamic 



becomes more complicated to understand for general values of h, since the "locality" property 
is lost. Figure [2] shows a situation when a maximum and a minimum disappear in one step 
(for this to happen, the area A has to be less than 2/i). Of course one can construct much 
more complicated examples. With this one, we can observe that the value of Uh inside [a;i,a;2] 
depends on the values of uq on both [a;i,a;2] and [x2,X3]. 







"jU) 




% 



Figure 2: Dynamics of a maximum step when h is not necessarily small. This figure shows the dynamic 
only inside the interval [xq , X4] . 



2.5 The dynamics of local step functions II. The continuous time case 



We now deduce the explicit evolution of strong solutions to the TVF introduced in Section 2.3 



for the Cauchy problem on M where the initial data coincide with a step function (see Remark 



2.7 below for the analysis of initial value problems with boundary conditions on intervals). The 
dynamic of step functions will then by used in the next section to deduce, by approximation, 
qualitative properties for general solutions. 



With the same notation as in Section 2.4.2, we consider the initial data given by 



N+l 



uo[x) 



y^ OkXh (^) inside /. 



(2.21) 



k=0 



Then, for h small and x £ I, we define 



u^{t,x) := (e+l--juih{x)+ f^-nju(^+i)ft(x) for any t G [eh,{e+l)h]. 



10 



Then it is immediate to check that, by (2.19), on the time interval [0,ti] with 



ti < min i la, — ai+i| min ||I,|, |/,+i|| > 
j=0,...,N U "^ -^ ' u Ji I JT- I J J 



the solution u is explicit, and independent of /i on /i U . . . U /at: 



N+l 



u^{t, x) = uq{x) + t ^ h/hXkix) on [0, ti] X /, 



fc=o 



with 



/3, 



klh 



0, 



2 

I4f 



2 

141 



if Ofc-i < Ofc < Ofc+i or if a^+i < a^ < a^.i 
if Q!fc > max{afe_i, afc+i} 
if Ofc < min{afc_i, afc+i} 



(2.22) 



for A; = 1, . . . , A^, and 



Po,eh 

PN+l,ih ■ 



>0, 
<0, 

> 0, 
<0, 



if ao < (^1 
if ao > «i 

if ttN > OAf+l 
if Oat < CKAf+l- 



By letting /i — )• 0, we find that i((x, t) remains a step function on /, its evolution is explicit on 
(/i U . . . U /at) , and on Iq and /at+i it is monotonically increasing/decreasing, depending on 
the value on /i and Ijy- 

This formula will then continue to hold until a maximum/minimum disappear: suppose 
for instance that 02 > max{ai,a3}. Then, after a certain time t'^, the value of u{t) on I2 
becomes equal to max{ai,a3}. Then, we simply take u{t\) as initial data and we repeat the 
construction. 

After repeating this at most N times, all the maxima and minima inside / disappear, and 
u{t) is monotonically decreasing/increasing on /. 

For instance, if I = M and the initial data is a compactly supported step function, then 
u = after some finite time T (which we call extinction time). On the other hand, if uq is an 
increasing (resp. decreasing) step function, then it will remain constant in time. 

Remark 2.7 The analysis done up to now can be extended to the case of suitable initial value 
problems on intervals with boundary condition. For instance, the dynamic of the Dirichlet 
problem is analogous to the one described above for the Cauchy problem with compactly 
supported initial data; we leave the details to the interested reader. Let us consider next the 
Neumann problem on some closed interval [a, b] = IqU. . .U/at+i. The dynamics on /lU. . .UIn is 
known by our analysis (which, as we observed before, is "local"). To understand the dynamics 
on Iq and /at, we go back to the time discretized problem: the Euler-Lagrange equations in 
this case are still (2.6), but with the additional Neumann condition Zhijj) = Zh{a) = 0. It is 



11 



easy to check that this last condition allows to uniquely characterize the value of Uh inside /q 
and In+1- 

For example, if uq = X]fe=i '^kXh with ai < . . . < aAf+i (i.e. uq is monotonically increasing), 
then 

u{t) =uo + t[ -rj-jXlo - 77 \XN+1 

Vl-'ol l-'Af+ll 

(i.e. the value on Iq increases, while the one on In+i decreases). This holds true until a jump 
disappears, and then one simply repeat the construction. We leave the details for the general 
case to the interested reader. 

2.6 Some properties of solutions to the TVF 

In this section we prove some local "regularity" properties enjoyed by solutions of the TVF 
(see also |131 115|). The key fact behind these results is that the TVF is contractive in any L"? 
space with q G [l,oo], cf. for example [H]. This contractivity property is not so surprising, 
since it holds also for the p-Laplacian for any p > 1. We are going to show that most of the 
properties which holds in the case of step functions, can be easily extended to the general case. 
An example is the following: 

Theorem 2.8 (Local continuity) Assume that uq is continuous on some open interval I . 
Then also the corresponding solution u{t) is continuous on the same interval I and the oscil- 
lation is contractive, namely 



supu(t) — inf u(t) =: osc {u{t)) < osc [uq). 



Proof. By the contractivity of the TVF in L°° , we have 

||u(t) -f(t)||oo < 11^0 -I'olloo (2.23) 

for any two given solutions u{t),v{t) corresponding to initial data Uo,Vo, with uq — vq £ L°° . 
Since uq is continuous in /, we can find a family of functions Uq such that Uq = uq outside /, 



Uq are step functions inside /, and ||no — Uq\\oo < e. Then by (2.23) we get 



\\u{t) — u^{t)\\oo < ||no — UqWoo < £ for all t > 

where u^{t) is the solution to the TVF corresponding to Uq, which is still a step function inside 
/. We now observe that the elementary inequality 

I osc(n(t)) - osc(u^(t))| < 2 \\u{t) - u%t)\\Loo^) < 2e (2.24) 

holds. Moreover, by looking at the explicit formulas for the evolution of u^{t) inside I, cf. 



Section 2.5, it is immediate to check that osc/(n^(t)) is decreasing in time. Hence 
osc (n(t)) < osc(n'^(t)) +2e < osc(uq) +2s < osc(no) + 4e. 

We conclude letting e — ;• 0. n 
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Remark 2.9 The above theorem still holds if uq is not continuous on /: in that case one has 
to replace sup and inf by esssup and essinf, and to prove the result one can use the comparison 
principle: if u'^{t) and u~{t) are the solution starting respectively from 

^ W '■- \ esssup -uo ifxGl; "" i^) '■- S essinf lio ii x e I; 

then u~{t) < u{t) < u^(t), 'U^(t) and u~{t) are both constant on /, and ||u^(t, x)—u^{t, x)\\icx,rj\ 
is decreasing in time. However, since we will never use this fact, we leave the details to the 
interested reader. 

2.7 Further properties 



Arguing by approximation as done in Theorem 2.8 above (using either the stability in L°° 
or simply the stability in L^, depending on the situation), we can easily deduce other local 
properties of the TVF, valid on any subinterval I where the solution u{t) is considered (we 
leave the details of the proof to the interested reader) : 

(i) The set of discontinuity points of u{t) is contained in the set of discontinuity points of no, 
i.e. "the TVF does not create new discontinuities" . 

(ii) The number of maxima and minima decreases in time. 

(iii) If no is monotone on an interval /, then n(i) has the same monotonicity as no on /. 
Moreover, if no is monotone on the whole M, then it is a stationary solution to the Cauchy 
problem for the TVF. 



(iv) As a direct consequence of Theorem 2.8 C'^'^-regularity is preserved along the flow for any 
a £ (0, 1]. (Similar results have been obtained for the denoising problem and for the Neumann 
problem for the TVF in [13].) Moreover, if no G W^'^{R), then n(t) e W^'^{R) (this is a 
consequence of the fact that the oscillation does not increase on any subinterval) . 

(v) If no G BViociR), ^ priori we do not have a well-defined semigroup. However, in this 
case no is locally bounded and the set of its discontinuity points is countable (see [H Section 
3.2]), and so in particular has Lebesgue measure zero. Hence, by classical theorems on the 
Riemann integrability of functions, we can find two sequences of step functions such that 
Uq < no < Uq and ||nQ — Uq ||i < e (the number of steps will in general be infinite, but 
finite on any bounded interval). Then, by approximation we can still define a dynamics, which 
will still be contractive in any L^ space. 

Behaviour near maxima and minima. We conclude this section by giving an informal 
description of the evolution of a general solution (excluding "pathological" cases) . 

Assume that no has a local maximum at xq. Then, at least for short time, the solution is 
explicitly given near xq by 

u{t,x) = min{no(x), /i(t)}, 

where the constant value h{t) is implicitly defined by 

[no(x) -h{t)] dx = 2t, 
'h 

Iq being the connected component of {no > h{t)} containing xq, see Figure |4J For a minimum 
point the argument is analogous. The dynamics goes on in this way until a local minimum 
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"merges" with a local inaxinium, and then one can simply start again the above description 
starting from the new configuration. 





Figure 3: Dynamic of TVF at maximum and minimum points. 



2.8 Rescaled flow and stationary solutions 

Let uq G BV{M.) n L-^(]R) be a nonnegative compactly supported initial datum. First we show 
that Uf) extinguishes in finite time, and we calculate the explicit extinction time. (Note that, 
even in general dimension, estimates from above and from below on the extinction time were 
already known, see for example [HEIIIH].) 

Proposition 2.10 (Loss of mass and extinction time) Letu{t) he the solution to the Cauchy 
problem in M for the TVF, starting from a non-negative compactly supported initial datum 
uq G L"'^(M). Then the following estimates hold: 



/ u{t, x)dx= / uo{x) dx-2t = 2{T - t) for all t > , 

Jr Jr 



and the extinction time for u is given by 

T = T{uo) = l [ 

^ JR. 



uq{x) dx . 



(2.25) 



(2.26) 



Proof. Arguing by approximation and using the stability of the TVF in L^, it suffices to 
consider the case when uq is a nonnegative step function. Assume that supp(uo) C [a, b] and 
that Uq jumps both at a and at b. Then, by the explicit formula for u{t), we immediately deduce 
that u{t) jumps both at a and b, with z{t,a) = 1 and z(t,b) = —1 (since u{t) is nonnegative 
as well). Hence 



d 
dt 



u{t, x) dx 



dxz{x, t) dx = z{b) — z{a) 



from which ( 2.25 )-( 2.26) follow. Q 



Remark. Let us point out that there is no general explicit formula for the extinction time 
when Uq changes sign. 
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The resettled flow. We now are interested in describing the behavior of the solution near the 
extinction time. To this end we need to perform a logarithmic time rescaling, which maps the 
interval [0,T) into [0,+oo), where T is the extinction time corresponding to the initial datum 
uq. We define 

w{s,x) = Y--^u{t,x) , Z{s,x) = z{t,x), s = riog f ^^^— ^ j , t = T i^l - e'"/^ 

(2.27) 
where u{t) is a solution to the TVF. Then 

w 
dsw{s,x) = d^Z+-, Z-D^w = \D^w\, w{0,x) = uo{x). (2.28) 

We observe that stationary solutions S{x) for the rescaled equation for w correspond to sepa- 
ration of variable solutions in the original variable, namely 

S T — t 

—dxZ = — provides the separate variable solution Uxit^x) := — — — S{x) . 

We need now to characterize the stationary solutions. To this aim, we first have to define the 
"extended support" of a function / as the smallest interval that includes the support of /: 

supp*(/) = inf {[a,5] | supp(/) C [a,b]} . 
Theorem 2.11 (Stationary solutions) All compactly supported solutions of the equation 

-a,Z = |, Z-DxS = \DxSl (2.29) 

are of the form 

2T 



S{x) = - XM](x), (2.30) 

— a ^ ^ 



with [a, h] C R. 



Proof. Let us assume that supp* 5 = [a, 6]. Since S is nonnegative we have Z[a) = 1, 
Z{h) = -1. We claim that -1 < Z < 1 on (a, b). 

Indeed, assume by contradiction that Z{xq) = 1 for some point xq G (a, h) (the case resp. 
Z{xq) = — 1 is completely analogous). Then, using again that S is nonnegative we obtain 
dxZ = — y < 0, which implies that Z = 1 on [a,xo]. Hence S = —TdxZ = on [a, xo], which 
contradicts the definition of supp*(5). 

Thanks to the claim, since Z ■ D^S = \DxS\ we easily deduce that \DxS\{a, b) = 0, that is, S 
is constant inside (o, b). To find the value of such a constant, we simply integrate the equation 
over [a,b], and we get 



/•b q pb 

/ dx = - dxZ{x)dx = Z{a)-Z{b) = 2. 

J a -'- J a 



and (2.30) follows. Q 
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Corollary 2.12 (Separate variable solutions) All compactly supported solutions of the TVF 
obtained by separation of variables are of the form 



T — t 

UTit,x) = 2—^X{a,b)ix) 



(2.31) 



where T > and [a, b] C R. 



Proposition 2.13 (Mass conservation for rescaled solutions) Let u{t) be the solution 
to the TVF corresponding to a nonnegative initial datum uq G BV(M.) H L^(M). Let w{s) be 



the corresponding rescaled solution, as in (2.27), then we have that 

w{s,x)dx = / uo{x)dx. 



(2.32) 



Proof. From (2.25), (2.27), and the fact that the extinction time is given by 2T = Jj^uodx, we 
deduce that 



w{s,x) dx 



T-t 



uq{x) dx — 2t 



.s/T 



/ no(x) dx - 2r (l - e"'*/^) = / uq{x) dx 
Jr ^ ^J Jr 



a 



Proposition 2.14 (Stationary solutions are asymptotic profiles) Letw{s,x) be a solu- 
tion to the rescaled TVF corresponding to a non-negative initial datum uq G BV{M.) n L-'^(M). 
Then there exists a subsequence s^ — )• oo such that w{sn, •) — > S in L^{I) as n —)• oo where S 



is a stationary solution as in (2.30). Equivalently we have that there exists a sequence of times 
T as n ^ oo such that 

u{tn,-) S 



t 



T-tr. 



T 



-^0. 



Li 



where S is a stationary solution. 



Proof. This is a weh known result, see e.g. Theorem 4.3 of [3] or Theorem 3 of [11 for the 
homogeneous Dirichlet problem on bounded domains. See also the book [5] . n 

Remark. From the above result we cannot directly deduce the correct extinction profile 



for the TVF, since there is not uniqueness of the stationary state, as Theorem 2.11 shows. 
Indeed a priori there may exists different subsequences such that the solution w approaches 
two different stationary states along the two subsequences. We shall prove in the next section 
that such phenomenon does not occur. 



2.9 Asymptotics of the TVF 

Here we want to characterize the asymptotic (or extinction) profile for solutions to the TVF 
in function of the non-negative initial datum uo € BV{M.) Cl L^(]R). 

Theorem 2.15 (Extinction profile for solutions to the TVF) Let u{t,x) be a solution 
to the TVF corresponding to a non negative initial datum uq G BV{M.) with supp*(iio) = [a, b], 
and set 



2 



^0(2;) dx . 
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Then supp(ti(i)) = [a, b] for all t G (0, T) and 



u{t, 



,X[a,b] 



Li(M) 



t^T 



-^0. 



(2.33) 



Remarks, (i) The above theorem shows to important facts: firstly, the support of the solution 
becomes instantaneously the "extended support" of the initial datum, which is the support of 
the extinction profile. Secondly, on [a, b] = supp*(iio) we consider the quotient u{t, x)/UT{t, x), 
where Ut is the separate variable solution Urit, x) = (T — t)S{x), and S{x) = 2-^j3^ (see ( 2.30| ) 
and ( 2.31| ) ) . It is interesting to point out that Ut is explicitly characterized in function of the 
extinction time (i.e 2 f uq) and of extended support of the initial datum. Then (2.33) can be 

rewritten as 

u{t,-) 



Uritr) 



Li(M) 



t^T 



-^0, 



(In literature this result is usually called convergence in relative error.) Equivalently, L^-norm 
of the difference decays at least with the rate 

||n(t,-)-t/T(t,-)||Li(M)<o(r-i)- 
In the next paragraph we will prove that the o(l) appearing in the above rate cannot be 



quantified/improved, so that the convergence result of Theorem 2.15 is sharp. 



(ii) The result of the above theorem can be restated in terms of the rescaled flow of Subsection 



2.8: \iw{s, x) is the rescaled solution corresponding to tt;(0, x) = uq (see 2.27), then the relative 



error 'w{s,x)/S — )• 1 as s — )• 00 in L-'^(supp*(uo))- 

Proof. The proof is divided into several steps. 

• Step 1. Compactness estimates. Let z{x, t) be associated to the u{t, x) (as in the definition 



of strong solution, see Section 2.3). We claim that 



\dxz{x,t) 



< 



2\\uo\ 



(2.34) 



The proof of this fact is quite standard in semigroup theory, once one observes that dtu{t, x) 
dxZu{x,t), Indeed, the homogeneity of the semigroup implies that 



\\dtuit,x)\\p< 



2||'Uo||p 



\/p G [1,00] . 



(2.35) 



Although the latter estimate is classical and due to Benilan and Crandall |12j , we briefly recall 
here the proof for convenience of the reader. Since u{t, x) is a solution to the TVF starting 
from u{0, •) = Mo(') ) it is then clear that ux{t, x) = Xu{X~^ t, x) is again a solution to the TVF 
starting from ux{0, •) = Amo(") , for all A > 0. Then 



u{t + h,x) — u{t, x) 



t + h 



ux{t,x) — u{t,x) = X ux{t,x) — u{t,x) 



(A ^ -l)uxit,x) + {ux{t,x) -u{t,x)) 
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where we have defined A := t/{t + /i) > 0. We now use the contraction property of the TVF 
in any L^— space to conclude that 

\\u{t + h,x) -u{t,x)\\p < (A^^ - l)\\ux{t,x)\\p + \\ux{t,x) -u{t,x)\\p 

< |A~^ - 1|||mo|Ip + \\ux{0,x) - u{0,x)\\p 

= |A"^ - l|||uo||p + |A- l|||no||p 



h h 

i TTh 



ml 



< 



2h. 



\uo\\p- 



Letting h -^ 0+, (2.35) follows. 



• Step 2. Stability up to T^ . Since uq £ BV{M), as in Subsection 2.7, Property (v), we 
can find two sequences of step functions Uq" < uq < u^ such that ||mq — Uq^||i < e. In 
particular, by the formula for the extinction time, we deduce that 

T-e/2< T{ul'') < r(no'+) < T + e/2 . 

Moreover, up to replacing Uq with minjuQ , ||no||ooXsupp*(Mo)}' ^^ have that 

supp*(tiQ ~), supp*(mq ) — )• supp*(mo) as e — )• 0. 

Since the evolution of step function is explicit, it is immediately checked that 

supp(n^'~(t)) = supp*(uq ~), supp(u'^'^(t)) = supp*(mq ) 

for t G (0, T — e/2] (indeed, for any subinterval I CC supp*(uQ~) where Uq'~ vanishes, u^'~{t) 
becomes instantaneously positive, see Figure [3J By the parabolic maximum principle, this 
implies 

supp*(uo") C supp(u(t)) C supp*(no'"'') for t € (0,T - e/2], 

which by the arbitrariness of e implies that 



supp('u(t)) = [a,b] Vt G (0,r). 



(2.36) 



• Step 3. Convergence of z and shape ofu{t) before the extinction time. By (2.36) and the fact 
that u{t) is nonnegative, we deduce that z{t, a) = 1 and z{t, b) = —1 for all t G (0, T). Estimates 
(2.34) imply that z{t,-) is uniformly boun ded in C^'^([a, 6]), hence it is compact in C^{[a,b]). 

Moreover, by (2.28) and Proposition 2.14 there exists a sequence Sk = Tlog i rp_^ j — )• oo (i.e. 

0. Hence, up to extracting a subsequence, we deduce that 



Efc — ^ 0^) such that dsw{sk,x 

lim z{T -ek) = l- [ 



A;— >c 



dy uniformly on [a,b], 



where S = -^^^ X[a,i3](.x) is a stationary solution, with a < a < j3 <b (observe that, by Step 2, 
the support can only shrink). 

• Step 4. Shape of u{t) before the extinction time. By Step 3 we know that z{T — e^) 
converges uniformly on [a, b] to the function 



zsix) 




ii a < x < a 
if a < X < /3 
if /3 < X < 6. 
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Hence, for k sufficiently large, there exist a < Oe,. < Pet < /? such that — 1 < 2;(T — e^) < 1 on 
[ae^,/3sj, z{T - Ek) > -1 on [a,Qej^], z{T - e^) < 1 on [/3ej^,6], and \a - ae^\ + |/3 - Pek\ -^ 
as e — 7- 0. Since z{T — e^) • Dxu{T — e^) = \Dxu{T — ek)\-, we easily deduce that u{T — e) is 
increasing on [0,0^^.], constant on [a^^, ,/3e^] and decreasing on [/3£j.,6]. 

• Step 5. Solutions with only one maximum point. Fix k large enough, and consider the 
evolution of u{t) on the time interval [T — efc,T]. By Step 4 and the discussion at the end of 



Subsection 2.7 the evolution of u{t) is explicit: 



u{t,x) = min{u(r — efc,x), h{t)} Vt G [T — efc,T) 
where h{t) > is implicitly defined by 

[^(r - Ek, x) - h{t)] dx = 2t. 



(2.37) 



Since supp(u(t)) = [a,b] (see Step 2), the above formula shows that the set [a(t),b{t)] where 



u{t) is constant expands in time and converges to [a, b]. Moreover, from (2.37) we easily obtain 

the estimate 

2 ; . X 1 

-7 <h{t)<— 

— a b — a 

for t close to T. Hence, the equation dtu = dxZ implies that 



< dxz{t) < 



b — a 



b — a 



on [a{t),b{t)]. 



(2.38) 



Since [a(t), 6(i)] — t- [a, 6], the uniform convergence of z{t) to zg is compatible with (2.38) if and 
only if 

( ■, 0, + b 2 
zs\x) = X on [a,oJ, 



i.e. the unique possible limiting profile is 2 T J^'''* , as desired. [] 



b — a b — a 

X[a,b] 

b-a ' 



2.10 Rates of convergence 

After proving asymptotic convergence to a stationary state, the next natural question is 
whether there exists a universal rate of convergence to it. As the next theorem shows, the 
answer is negative. 

Before stating the result, let us make precise what do we mean by decay rate. 

Definition 2.16 Let ^ : [0,oo) — t- [0, 00) be a continuous increasing function, with ^(0) = 0. 
We say that ^ is a rate function if, for any solution u{t) of the TVF, 



u{t) S 
T -t ~ T 



LHI) 



< ^{T — t) for any t close to the extinction time T . 



(2.39) 



The following result shows that there cannot be a universal rate of convergence to any station- 
ary profile. 
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Theorem 2.17 (Absence of universal convergence rates) For any rate function S, ■ [0,oo) 
[0,00) , there exists an initial datum uq G BV(R), with supp*('Uo) = [0, 1], such that 



2i{T -t) < 



u{t) 



T-t 



2X[o,i] 



for any <T — t < 1. 



LHi) 



(2.40) 



Proof. Let us fix a rate function ^ : M — )■ M. It is not restrictive to assume that ^ is strictly 
increasing, and tliat ^(s) > s for any s £ [0, 1]. 

Let ^~^ = [0,00) — ;■ [0,00) be the inverse of ^, so that ^^^{s) < s for any s £ [0,1], and 
choose the initial datum uq to be 



uo{x) 



1, 
co^^l-x), 



if < X < ^ 
if I < X < f 
if f < X < 1 



(2.41) 



with Co := l/C"Hl/4) (so that uo{l/A) = uo(3/4) = 1). By Theorem [2J5| we know that the 
solution u{t) corresponding to uq extinguish at time 

1 1 /"^ 1 

4-^=270 ^°(^)^^^2 

and that u{t)/{T — t) converges strongly in i^([0, 1]) to S/T = 2x[o,i] as t — )• T. First we prove 
that the L°°-norm satisfies the bound 



2{T-t)< ||n(t)||oo<4(r-t). 



(2.42) 



The first equality follows by the loss of mass formula (2.25) and Holder inequality on [0, 1]: 

2{T-t) = \\u{t)\\i< \\u{t)\\^. 
The second inequality follows as we know the explicit behaviour of the solution around maxi- 



mum points (see the end of Subsection 2.7), namely 

n(t, x) = min{uo(x), h{t)} 

where h{t) > satisfies 

/ [uo(x)-/i(t)] dx = 2t. (2.43) 

Then u{t) is constant on an interval of the form [Q;(t), /3(t)], with a{t) — )• 0"*" and /3(t) — )• 1~, 
see Figure |4j 

Hence 



/■I /•3/4 

2{T - t) = \\u{t)\\i = j u{t,x)dx> I u{t,x)(ix 

Jo Jl/2 



\U 



and so \\u 



< 4(T — t) , and inequality (2.42) is proved. 



Now, let ao(*) G [0)Ct(0] ^^'^ /3o(0 ^ [/^(O)!] be the unique points such that ti(t, ao(i))/(T — 
t) = u{t, /3o(t))/(T — t) = 2 (such points exists thanks to the lower bound ||M(t)||oo > 2{T — t), 
see also Figure |4J). Let T^q (resp. T^^) be the rectangular triangle with height {0} x [0,2] 
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Figure 4: Left: Dynamic of u{t): black: uo(x) , blue: u{t^x) ^ red: u{t -\- h) . Right: Reseated 
dynamic: black: U{){x) (dashdot) and u^/T (cont.), blue: S{x) = 2x[o,i], red: u{t^x)/{T — t)^ green: 
definition of q;o(^)i /5o(0 ^^*^ ^ao ■> ^(i^ • 

(resp. {1} X [0,2]) and basis [0,a(t)] x {2} (resp. [/3(t),l] x {2}), as depicted in Figure [4| 
Denote by IT^qI, iT^ol their measure. Then, since ^~^(s) > 5 we easily obtain the estimate 



u(t) 


S 




_ r 


u{t,x) 


T -t T 


Li([0,l]) -^0 


T-t 


Jo 


T 


-t 


dx + 


Ja{t) 



-2x[o,i](2;) 
u(t, x) 



> 



y{t) 







u{t,x) 



T-t 



dx + 



m 



T-t 

u{t,x) 



dx 



dx + 



m 



u{t,x) 



T-t 



dx 



T-t 



dx 



> J"'" (^2-'f^] 4.^ " 



Mt) 



T-t 



>|Taol + |T/^ 



2ao(t). 



To estimate ao{t) from below, we observe that on [0, a{t)] we have that u{t, ao^t)) = uo{ao{t)) 
coC~H"o(0), so 

ao{t)=c(-{T-t) 
\co 

Now, recalhng that ^ is strictly increasing and .^(1) > 1, we get ^(2) > 1/4, or equivalently 
2/co = 2/^-1(1/4) > 1. Hence ao{t) > ^(T - t), which concludes the proof. D 

Remark. The above Theorem shows that there cannot be universal rates of convergence. A 
similar construction will provide (nontrivial) initial data for which the convergence is as fast 
as desired. 



Theorem 2.18 (Fast decaying initial data) For any rate function ^ : [0,oo) — ;• [0,oo), 
there exists an initial datum uq G L^{I) such that the corresponding solution u(t) satisfies 



T-t 



2X[o,i] 



Li(/) 



<^(8(r-t)), for anyO<T-t<l 



(2.44) 
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Proof. Fix a rate function .^ : M — )• M, which is continuous, increasing, ^(0) = 0, and ^(s) < s. 

Let ^~^ = [0,00) — )■ [0,00) denote its inverse, and define uo as in (2.41), see Figure [sj Then 
an analysis analogous to the one done in the previous Theorem proves the result. We leave 
the details to the interested reader. □ 



0,2- r 
I 
i 



0,6 o.s 





Figure 5: Left: Dynamic of u{t): black: uo(a;),blue: u(i,a;),red: u{t + h,x) . Right: Reseated 
dynamic: black: UQ{x),hhxe: S(x) = 2x[o,i], red: u{t,x)/{T — t), green: definition of Q'o(i),/3o(0 ^i^d 

Tan ' ^I3a ■ 



3 Solutions to the SFDE and solutions to the TVF 

As explained in the introduction, TVF and SFDE are formally related by the fact that "u 
solves the TVF if and only if D^u solves the SFDE". 

In order to make this rigorous, we need first to explain what do we mean by a solution of 
the SFDE, and then we will prove the above relation by approximating the TVF with the 
j)-Laplacian and the SFDE by the porous medium equation. 

The notion of solution we consider for the SFDE is the one of mild solution. More precisely, 
since the multivalued graph of the function r 1— t- sign(r) is maximal monotone, by the results 
of Benilan and Crandall [TT], there exists a continuous semigroup Sq : L^(M) — t- L^(M) such 
that S^vq G C([0, 00); -L^(M)) is a mild solution of the SFDE. To be more precise, let if he a 
maximal monotone graph in M (see [lOj ) and consider the problem 



u{0,x) = uo{x) , 



in P' ((0,00) X 

xGM 



where the first equation is meant in the sense that 



Ut = Aw 



m 



V'{{0, 



00 X 



with 'w{t,x) S ip{u{t,x)) 



a.e. t, X G 



(3.1) 



(3.2) 



We now recall the celebrated results of Benilan and Crandall |llj adapted to our setting, 
namely it is sufficient to consider (p{r) = sign(r). 

Theorem 3.1 (Benilan-Crandall, |llj ) Under the running assumptions, the following re- 
sults hold true: 
(i) There exists a unique solution u G C([0, 00); L"^(M)) n L°°([0, 00) x M) corresponding to the 
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initial datum uq G L-'^(M) nL°°(]R) such that (|3.1[) and (3.2) are satisfied. 



(ii) Let Un G C{[0,oo); L^(M.)) n L°°([0,oo) x M) 6e solutions of (3.2) corresponding to the 
sequence c/?„ : M — )■ M, n = 1, 2, . . . o/ maximal monotone graphs in M. Assume that G V'n(O), 



lim V'nl'") = v{t) for r G 



and 



lim ||no,n --uoIIlm 



0. 



(3.3) 



Then Un ^ u in C([0,oo); L-'^(M)') , where u is the solutions of (3.2) corresponding to ip 



Now, let 5™ be the semigroup associated to the FDE equation 

dtv = ^{v"^). 



3.1 



^yn ij- converge to the graph of r i— t- sign(r), 
|r|"~ r to guarantee that we have 



le function r i— t- r"* 

with the simple choice (/9„(r 

0+) of SY'VQ to S'^vq in C([0, oo); L1(M)), for any initial datum 



Since the graphs of t 
we can use Theorem 
convergence (indeed as m 
VQ G L\ 

On the other hand, we can consider the p-Laplacian semigroup T[ for p = 1 + m. It is well 
known (see for example HQj chap. 4, or [5] chap. 5) that if uq G W^'P(R), then T[uo G W^'^iR), 
so that as p — )■ 1"*", strong solutions to the p-Laplacian converge to strong solutions to the TVF. 
(A detailed proof of this fact in dimension n > 1 can be found for instance in [Ij, pg. 138-142, 
in the framework of the Dirichlet problems on bounded domains, but it can be easily adapted 
to other problems, including the Cauchy one.) Hence, if uq is a smooth compactly supported 
function, Tfuo — )• T^uq in C([0, cxd); L^(M)) as p — )• l"*", where T/ denotes the TVF-semigroup. 

Moreover, if p = 1 + m, we have that dx(TfuQJ solves (in the distributional and semigroup 
sense) the FDE with initial datum dxUo, i.e. dx (Tf mq) = S]^ {dxUo) ■ Hence, by letting ?n, — )• 0+, 
we recover such a relation in the limit p = 1 and ?n, = by what we said above. We can 
summarize this discussion in the following: 

Theorem 3.2 Assume uq is a smooth compactly supported function. Let 1 < p < 2, and 
m = p — 1. Then the following diagram is commutative: 



Tf-uo G W^^P{1 



p 



r 



T}uq G P7^'^(I 



dx 



sr{dxUo)eL'-^^iI 



S?{dxUo)eLH 



m ^ 0^ 
Note that the convergence in meant in the sense of distributions 

At this point it is worth noticing that the vector field z associated to the solution u of the 



TVF as in Section 2.2 and the function w associated to the solutions to the SFDE as in (3.2) 



are the same (just by letting vq = OxUq and v = dxu) 

dtu = dxZ > dfV = dtdxU = dxxZ = dxxW . 



Measures as initial data. Once the correspondence between TVF and SFDE is established 
for smooth initial data, by stability in L^ of both semigroups it immediately extends to uq G 
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Ty"'^'^(M), and then by approximation to BV{M.) n L^(R) initial data. However, at the level of 
the SFDE this would correspond to finite measures vq such that J^ vo{dy) G -L"^(M), which is 
possible if and only if f_^ vo{dy) = 0. Actually, this class of data correspond exactly to the 
one for which there is extinction in finite time (as this is the case for L^ initial data to the 
TVF). 



To remove this unnatural constraint on vq, we observe that, by Subsection 2.7 Property (v), 
the TVF defines a contractive semigroup also on initial data which are only in BViocO^)- ^^ 
particular, the TVF is well-defined on data of the form ^0(2;) = J_^vo{dy), where vq is a 
(locally) finite measure on M. Hence, this allows to define measure valued solutions of the 
SFDE as dxTf{uo), and this notion coincides with the one that one would get by considering 
weak* limit of L^ solutions. 

Summing up, we have shown that: 



If Vq £ L^(M), the unique mild solution of the SFDE of Theorem 3.1 is given by 

Sfvo = d.(^n(^J' voidy))y (3.4) 



Using (3.4) we can uniquely extend the generator S^ to measure initial data (actually, since 



the semigroup T/ is well-defined on L^(]R), one could even extend the SFDE to distributional 
initial data in W-^^'^{R)). 

3.1 The 1-dimensional SFDE 



By (3.4), the dynamics of the SFDE can be completely recovered from the one of the TVF. 
We begin by illustrating some basic properties, and then, instead of trying to give a complete 
description of the evolution (which, by the analysis of the TVF done in the previous sections, 
would just be a tedious exercise), we prefer to briefly illustrate the evolution of solutions of the 
SFDE in some simple but representative situations. As pointed out in the previous section, 
we are allowed to consider measures as initial data for the SFDE, keeping in mind that the 
distributional j;-derivative of a solution to the TVF is a solution to the SFDE. 

In the same way as step functions allowed us to understand the dynamics of the TVF, we 
start by considering the dynamics of sum of delta masses, which is directly deduced from the 
one of step functions for the TVF. 

Example 1. Delta masses as initial data. Let us assume that uq = ^i=iO'i^Xii with 
xi < • • • < xn- Then, for i > small (the smallness depending on the size of \ai\) 

N 



v{t) = y^^ai{t)5x 



i=l 




if sign(ai_i) = sign(ai+i) = sign(aj) 
sign(ai)(|ai| - At) ^ , if sign(ai_i) = sign(ai+i) = - sign(ai) (3.5) 

sign(ai)(|ai| - 2t) , if sign(aj_i) sign(aj+i) = -1, 
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Figure 6: black: vo{x) , blue: v{t,x), t < to, red: v{ti,x) is a stationary state. 



where we use the convention sign(ao) := sign(ai) and sign(aAr+i) := sign(aAr). This formula 
holds true until one mass disappear at some time t'^ > 0, and then it suffices to v{t'i) as initial 



data and repeat the construction (compare with Subsection 2.5). 

We observe that if all Oj have the same sign, then vq produces a stationary solution v{t, x) = 
vo{x). Moreover, the total mass is conserved under the dynamics. In particular, v{t) extin- 
guishes in finite time if and only it has zero mean, i.e. ^^ Oj = (however, there is no simple 
formula for the extinction time, see the remark after Proposition |2.10|. 



General properties. Arguing by approximation (or again using the direct relation with the 
TVF), as a consequence we have the following properties of the SFDE flow: 

(i) Nonnegative initial data. Let wq > be a locally finite measure, and define ^0(3;) := 
/^oo^o^^^) — ^' Since uq is monotone non-decreasing, it does not evolve under the TVF, cf. 
(iii) in Section 2.7 Hence vq is a stationary solution to the SFDE. (Actually, since monotone 
profiles are the only stationary state for the TVF, the only stationary solutions for the SFDE 
are nonnegative/nonpositive initial data.) 

(ii) Only zero mean valued initial data extinguish in finite time. If vq is a finite measure, v{t) 
converges in finite time to a stationary solution v such that J^v{dy) = f^vo{dy). Moreover, 
V = (i.e. Vq extinguish in finite time) if and only if J^ vo{dy) = 0. 



Example 2. Interaction between a delta and a continuous part. 

Let Vq = vo+a6o where vq G C(M) is positive exception made for an interval [0, xq] as depicted 
in Figure^ The zero set of vq is Z{vq) = [0, a] U {xq}, and a Dirac's delta with mass a is put 
at X = 0. We assume a is much smaller than J^" \vo\. The evolution basically changes at two 
steps: 

• < t < to- The delta mass starts to lose its mass a by a factor 2t until time to = ck/2 
(when it extinguish). This mass is compensated by a "gain of mass" of vq: its zero set near 
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a starts to move to the right, and at time t it is at a position zi{t) characterized by the 
fact that J" ^^ It'ol = 2t (cf. the blue area 2t in the Figure 6). On the other hand, the 
isolated point xq G Z{vq) starts to "expand", creating a zero set[z2{t) , Z3{t)] G Z{u{t)), with 
Iz^d) l^ol ~ /^^ l^ol ~ ^^' This expansion of the zero set corresponds to the creation of flat 
parts at the level TVF (see Figure [s]). 

• to < t <ti. At time to the delta disappears, and u{to,x) is a piecewise continuous function 
on M which is positive outside [0,2:3(^0)], it is zero on Z{u{tQ,x)) = [0, zi(to)] U [22(^0); -23 (to)]; 
and is negative on (zi(to), Z2{to)). Starting form this time, the free boundary expands on both 
components in the same way as described above, until some time ti > to when zi{ti) = Z2{ti). 
Observe that the loss of mass is such that Bq = Bi + B2 (see Figure pi). 

• Reaching the stationary state infinite time. Since u{ti) > 0, the solution becomes stationary 
and u{t) = u{ti) for t > ti. 

3.2 SFDE vs LFDE. 



Let us go back to the fast diffusion equation (1.1), and assume that vo ^ (so v{t) > for all 



t > 0). We remark that, changing the time scale t 1— t- mt, the above equation can be written 
in two different ways which lead to two different limiting equations: more precisely, setting 
p{t,x) = v{t/m,x), 

dtv = A[v"^) ¥ dtv = A[sign{v)) 

m— >• 0+ 

dtp = div(y9'"-iVp) > dtp = div (p-^Vp) = A(log(p)) . 

Observe that the evolution of p{t, x) on the time interval [0, T] corresponds to the evolution of 
v{t, x) on the larger time interval [0, T/m], and the diffusion of v is slower than the diffusion of 
/? by a factor 1/m. So, when analyzing the limit as ttt. — )■ 0^, one gets two different limits, and 
the evolution of p(t, x) on the time interval < t < T corresponds to the evolution of v{t, x) 
on the time interval < t < 00 for every T > 0. This means that the solution to the SFDE 
corresponds to an evolution "infinitely slower" than the solution to the LFDE. 

We saw that the Cauchy problem for the SFDE gives rise to a trivial dynamics on non- 
negative initial data. However, the problem becomes non-trivial if we consider for instance the 
Dirichlet problem for the SFDE on a closed interval / with zero boundary conditions. Indeed, 
by "integrating in space" such a solution, we obtain a solution to the TVF on I with Neumann 



boundary conditions, whose dynamics on step functions has been described in Remark 2.7 In 



particular, from the example given there, we can explicitly find the dynamics of a finite sum of 
positive deltas: if vq = Yli=i o-i^xi with Oj > and 2:1 < ... < xjy, then v{t) = vq — t[6x-i +^xn] 
until one delta disappears, and then one simply restarts from there. By approximation, we 
see that positive initial data extinguish in finite time, and the extinction time is given by 
2/R^o(dy). 

This fact and the above discussion suggest that solutions to the LFDE should extinguish 
instantaneously, i.e. p(t) = for any t > 0. As one can deduce by the results of Rodriguez and 
Vazquez [22j, this is actually the case. (See also [HI EU |22] and the book [5B] for a complete 
theory of the logarithmic diffusion equation in one space dimension for positive initial data.) So 
the above discussion is correct, and gives an heuristic explanation for this immediate extinction 
phenomenon. 
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